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Soybean dwarf virus (SbDV) causes serious dwarfing, yellowing and sterility in soybean (Glycine max). The 
soybean cv. Adams is tolerant to SbDV infection in the field and exhibits antibiosis to foxglove aphid 
(Aulacorthum solani), which transmits SbDV. This antibiosis (termed "aphid resistance") is required for tol- 
erance to SbDV in the field in segregated progenies of Adams. A major quantitative trait locus, Rasol, is 
reported for foxglove aphid resistance. Our objectives were to fine map Rasol and to reveal whether Rasol 
alone is sufficient to confer both aphid resistance and SbDV tolerance. We introduced Rasol into cv. 
Toyomusume by backcrossing and investigated the degree of aphid antibiosis to foxglove aphid and the de- 
gree of tolerance to SbDV in the field. All .KasoZ-introduced backcross lines showed aphid resistance. Inter- 
estingly, only one .Kasoi-introduced backcross line (TM-1386) showed tolerance to SbDV in the field. The 
results demonstrated Rasol alone is sufficient to confer aphid resistance but insufficient for SbDV tolerance. 
Tolerance to SbDV was indicated to require additional gene(s) to Rasol. Additionally, Rasol was mapped 
to a 63-kb interval on chromosome 3 of the Williams 82 sequence assembly (Glymal). This interval includes 
a nucleotide-binding site-leucine-rich repeat encoding gene and two other genes in the Williams 82 soybean 
genome sequence. 
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Introduction 

Soybean dwarf virus (SbDV), a member of the Luteoviridae, 
is known to infect to soybean (Glycine max), red and white 
clovers (Trifolium pratense and T. repens) and other legume 
species (Tamada 1970). Infection of soybean by SbDV is 
widespread in northern Japan (Tamada 1975) and the United 
States (Thekkeveetil et al. 2007). In nature, SbDV is persis- 
tently transmitted by several aphid species (Gildow et al. 
2000). In northern Japan, the foxglove aphid (Aulacorthum 
solani) transmits SbDV (Mikoshiba et al. 1995) from 
SbDV-infected red and white clovers to soybean during 
spring to summer. Infection by SbDV causes serious dwarf- 
ing, yellowing and sterility, and leads to significant yield 
losses in soybean production. To avoid damage induced by 
SbDV, soybean breeders have attempted to introduced 
SbDV resistant/tolerant germplasm. An old American culti- 
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var, 'Adams' (PI 548502), was found to show lower SbDV 
symptom incidence in the field. Adams suffers SbDV infec- 
tion to a limited extent. However, the virus infection rate and 
symptom incidence are lower than other cultivars (Tanimura 
and Tamada 1976). In this paper, we termed this partial re- 
sistance to SbDV in the field "tolerance". Interestingly, the 
number of foxglove aphids observed on Adams is lower than 
that on other cultivars in the field (Kanehira 1978). More- 
over, young Adams plants exhibited antibiosis to the fox- 
glove aphid in a no-choice test (Jinno et al. 1997). We 
termed this aphid antibiosis "aphid resistance". Resistance 
to foxglove aphid was required for SbDV tolerance in the 
field when SbDV tolerance and aphid resistance were inves- 
tigated in segregated progenies of Adams (Araki et al. 

1997) . 

Aphid resistance is known among many plant species, 
including tomato (Solanum lycopersicum) (Rossi et al. 

1998) , wheat (Triticum aestivum) (Liu et al. 2005), barley 
(Hordeum vulgare) (Mittal et al. 2008), melon (Cucumis 
melo) (Brotman et al. 2002), barrel clover (Medicago 
truncatula) (Klingler et al. 2005), maize (Zea mays) (So et 
al. 2010) and soybean (Kim et al. 2010a, 2010b, Zhang et al. 
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2010). Some aphid resistance genes have been cloned or fine 
mapped. In most cases, the resistance genes or the candidate 
of the resistance gene are a typical R-gene containing a 
leucine-rich repeat (Dogimont et al. 2010, Kim et al. 2010a, 
2010b, Klingler et al. 2005, Milligan et al. 1998, Rossi et al. 
1998, Smith and Boyko 2007). Thus aphid resistance in 
most plants seems to be controlled by the same mechanism 
as for other pest and disease resistance. 

A previous report showed that aphid resistance in Adams 
is controlled by a quantitative trait locus (QTL) named 
Rasol (Kamiya et al. 2008), which has a relatively high con- 
tribution rate (r 2 = 0.32). Using 204 RILs, Rasol was detect- 
ed on the 3.7 cM region between Satt009 and Satt530 on 
chromosome 3 with a lod score of 16.0. The QTL analysis 
showed Rasol was essential for aphid resistance, but it is not 
clear whether Rasol is solely responsible for the aphid resis- 
tance. Similarly, based on previous reports, Rasol was de- 
duced to be required for the SbDV tolerance phenotype of 
Adams, but it is not clear whether Rasol alone is sufficient 
to confer SbDV tolerance. In this study, we fine-mapped 
Rasol and investigated whether Rasol alone is sufficient to 
confer both aphid resistance and SbDV tolerance phenotype 
in Adams. 

Materials and Methods 

Plant material 

Soybean cv. Adams (PI 548502) shows both foxglove 
aphid resistance and SbDV tolerance in the field (Jinno et al. 
1997, Kanehira et al. 1978, Tanimura and Tamada 1976). 
Shokukei-10 is a breeding line which shows both foxglove 
aphid resistance and SbDV tolerance in the field (Fig. 1). 
Shokukei-10 was derived by crossing Adams with Mechita, 
a Russian cultivar (Fig. 1A). We introduced the aphid re- 
sistance QTL Rasol by backcrossing to Toyomusume, 
which is a commercial cultivar in northern Japan. Four back- 
crosses were performed using Toyomusume as the recurrent 
parent and Shokukei-10 as the donor of Rasol. A total of 9 
backcross lines were developed from 9 BC 4 Fi plants. The 
Toyomusume backcrossed lines (TM-lines) were named 
TM-1382, TM-1386, TM-1387, TM-1388, TM-1394, TM- 
1397, TM-1402, TM-1404 and TM-1409. The simple se- 
quence repeats (SSR) markers Sattl52, Satt009 and Satt530 
were used to introduce Rasol based on a previous report 
(Kamiya et al. 2008). Primer sequences and positions of 
SSR markers on the genome sequences are described in 
SoyBase (www.soybase.org). We started the backcrossing 
by crossing Toyomusume to F, plants generated from a 
cross between Shokukei-10 and Toyomusume. Then, 
Toyomusume was crossed to a population consisting of a 
mixture of 4 independent BCiFi plants. To generate BC 3 F| 
and BC 4 Fiplants, 12 randomly selected individual BC 2 F| 
and BC 4 Fi plants were used. Thus, the final 9 BC 4 lines 
were partially independent. We used mapping populations 
C1905 and C1950, which were originally made for a breed- 
ing program and therefore have complex pedigrees. The 
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Fig. 1. Schematic representation of development of experimental 
lines. The plus (+) or minus (-) signs in parentheses represent the pres- 
ence or absence of QTL Rasol, aphid resistance, and SbDV tolerance, 
respectively, derived from Adams. Question marks represent "deter- 
mined in this paper". (A) TM-lines are Toyomusume-backcrossed 
lines with the putative aphid resistance QTL Rasol. (B) CI 905 and 
CI 950 are mapping populations. CI 905 and CI 950 were originally 
made for a breeding program, but we used them for fine mapping. 
C1950 was derived from a 3-way crossing of F! plants of C1905 and 
Toyoharuka. 

population CI 905 was derived by crossing the population 
derived from mixing the 9 TM-lines described above in the 
F 4 generation to an experimental line, Shokukou-05405-1-1 
(Fig. IB). Shokukou-05405-1-1 is a breeding line derived 
from backcrossing Tyokomachi, a cultivar in northern 
Japan, to an Indonesian cultivar, Wilis. The mapping popu- 
lation C1950 was generated by crossing Fj plants of C1905 
and Toyoharuka, a cultivar in northern Japan. Crossing 
Toyoharuka was for breeding purposes and not essential for 
mapping. Shokukou-05405-1-1 and Toyoharuka do not ex- 
hibit aphid resistance. 

Evaluation of aphid resistance 

For nine TM-lines (TM-1382, TM-1386, TM-1387, TM- 
1388, TM-1394, TM-1397, TM-1402, TM-1404 and TM- 
1409), a no-choice test was performed according to a previ- 
ous report (Jinno et al. 1997) to determine the degree of 
aphid antibiosis. Foxglove aphids (Acyrthosiphon solani) 
were collected atNaganuma, Japan (43.017°N, 141.695°E), 
and were maintained on fresh leaves of the aphid-susceptible 
soybean cultivar Chusei-hikarikuro. Three soybean seeds 
were sown in unglazed pots (height 12 cm; diameter 10 cm) 
filled with potting soil. Seedlings were thinned to one per 
pot to synchronize the plant growth stage. A single seven- 
day-old aphid was settled per plant when the first true leaf 
was 3-5 mm long (cotyledons were not completely opened). 
After settlement, plants were grown in a plant incubator 



620 Ohnishi, Miyake, Takeuchi, Kousaka, Hiura, Kanehira, Saito, Sayama, Higashi, Ishimoto, Tanaka and Fujita 



Table 1. Primer sequences and location of markers on chromosome 3 of the Williams 82 soybean genome sequence assembly (Glymal) 



Marker Name 


Location" 


Forward primer 


Reverse primer 


Gm03-01 


4,811,664 


AACCGTCGCTTAATCGTGATA 


CCCTGAAAGCATAAAAAGTGG 


Gm03-03 


5,448,257 


CATAAACTTTGGTCATGCTTAACTG 


ACCTTGACCAAACATTCTTGA 


Gm03-06 


4,683,586 


CTTGTATATATTACATTAAAATTGGAG 


TGATGAAGAGTGACTACAGTTGAG 


Gm03-07 


4,697,806 


AAAAACCTTGATTTGTTCTTCA 


CAAAGGTTGTTACAGAAATTCAA 


Gm03-09 


4,613,933 


CCATTCTAAGATGATTATTTACCTG 


GACGCATGAAGACGATTTT 


Gm03-10 


4,637,354 


CGGACAACCAAATCTCAAT 


GA A AGC ATGG AATA A AC GA A 


Gm03-ll 


4,661,084 


GTGGCCACATGAAAATACAG 


TCCTGGTGGACTTGCTAATA 


Gm03-12 


4,724,159 


AAAAAGGCTAAAAGTATATTTCTGT 


TGATACCATGATGTTTATATTTTT 



" Nucleotide numbers from the end of chromosome to the forward primer. 



maintained at 15°C with a 16 h/8 h (light/dark) photoperiod. 
Each plant was covered with a mesh cage. The number of 
aphids on each plant was counted 14 days after settlement. 
The same procedure was adopted for fine mapping of Rasol. 

Evaluation of SbDV tolerance in the field 

To evaluate SbDV tolerance, we selected 3 TM-lines 
(TM-1394, TM-1382 and TM-1386) that differ in degree of 
aphid antibiosis. Field evaluations of SbDV tolerance were 
performed in Date City, Japan (42.505°N, 140.889°E), a re- 
gion with a severe SbDV endemic infection. There are two 
virus strains SbDV-D and SbDV-Y, which cause yellowing 
and dwarfing symptoms in soybean (Tamada 1973). Both 
strains are observed in Date. In this region, no viral disease 
except SbDV was observed in soybean, and thus the dwarf- 
ing or yellowing symptoms were likely caused by SbDV. 
Soybean lines were sown on 13 May 2008 and 8 May 2009. 
Soybean fields were managed according to the local custom 
except that pesticide applications for control of aphids were 
omitted. Endemic SbDV symptoms on soybean were first 
observed about 50 days after sowing. Soybean plants with 
dwarfing or yellowing symptoms were counted and the inci- 
dence rate against the total number of the tested plants was 
calculated as "virus symptom incidence". Four replications 
of each line were grown and each 1.32 m 2 plot contained 22 
plants. Virus symptom severity was scored according to a 
previous report (Mikoshiba et al. 1990) on 3 November 
2009, 113 days after sowing. 

Mapping and marker assay 

From the mapping populations described above, a total 
of 669 F 2 plants for C1905 and 576 F 2 plants for C1950 
were screened with the SSR markers Satt009 and Satt530, 
which flank the QTL Rasol (Kamiya et al. 2008). We se- 
lected 59 plants for C1905 and 34 plants for C1950 in which 
a recombination event between the markers was confirmed. 
The harvested plants were used to produce the next genera- 
tion in a greenhouse and one homozygous plant per re- 
combination event was selected and grown to produce seed. 
For further genotyping between Satt009 and Satt530, SSR 
were identified from the Williams 82 soybean genome 
sequence GlymaOl (Schmutz et al. 2010) and primers were 
designed. After PCR amplification, the SSRs that showed 
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Fig. 2. Aphid resistance phenotype of ^iaso./ -introduced Toyomusume- 
backcrossed lines (represented by TM lines). (A), (B) and (C) show 
the results of three independent experiments. Shokukei-10 is an 
Adams-derived breeding line that shows both foxglove aphid-resistant 
and SbDV-tolerant phenotypes. Vertical axis indicates the number of 
aphids per plant 14 days after settlement in a no-choice test. TM lines 
highlighted in a box were evaluated for their virus tolerance phenotype 
(see Figs. 2, 3). Vertical bars indicate the standard error of the means 
(n = 8). Static comparisons were carried out independently for each ex- 
periment. Grouping by the square brackets below the line names indi- 
cates there was no significant difference within the group at the 5% 
level using the Tukey-Kramer multiple comparison test. 



polymorphism between Adams and Shokukou-05405-1-1 or 
Toyoharuka, the parents of the mapping populations, were 
used as SSR markers. The position of the SSR markers on 
chromosome 3 and the sequences of the primers are shown 
in Table 1. The genotype of Satt641 and the newly devel- 
oped SSR makers were assessed for 42 selected lines from 
93 plants in which the recombination event was confirmed. 
For genotype data, 1 5 informative lines were used to evalu- 
ate the degree of aphid antibiosis. For screening of the map- 
ping population, DNA extraction and SSR marker assay 
were performed as described previously (Sayama et al. 
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Fig. 3. SbDV tolerance phenotype of Toyomusume-backcrossed 
(TM) lines in the field. Virus symptom incidence in 2008 (A) and 2009 
(B) are shown. Dwarfing or yellowing symptoms induced by SbDV 
were examined in the field. The X-axis represents the assessment date. 
Vertical bars indicate the standard error of the means (n = 4). The stan- 
dard errors for TM-1382 and TM-1394 are of the same level as 
Toyomusume and its display are omitted. The p-values by Dunnett's 
test are indicated on the right side of each panel. Dunnett's test was 
performed for the final observations in each year using Toyomusume 
as the reference. 
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Fig. 4. Virus symptom severity of Toyomusume-backcrossed (TM) 
lines 118 days after sowing in 2009. Dwarfing or yellowing symptoms 
induced by SbDV in the field were assessed using the scoring system 
of Mikoshiba et al. (1990). Each plant was scored on a scale between 0 
(no symptom) and 4 (severe symptoms) with increments of 1. The 
symptom severity was calculated by [100(Sum of the score of plant)/ 
4(total number of plant)]. Vertical bars indicate the standard error of 
the means (« = 4). The p-values by Dunnett's test using Toyomusume 
as the reference are indicated at the top of the bars. 



Toyomusume 14 days after aphid settlement (Fig. 2). How- 
ever, the degree of aphid antibiosis varied among lines. 
TM-1386 and TM-1394 exhibited a similar number of 
aphids per plant as compared to Adams and Shokukei-10, 
from which the aphid resistance originated, whereas the 
other TM-lines exhibited a higher number of aphids than 
Adams did (Fig. 2). However, we concluded that the aphid- 
resistance QTL Rasol solely confers the aphid resistance 
phenotype when it is introduced by backcrossing. 



2011). To obtain the sequence and gene annotation data of 
Williams 82, Phytozome v7.0 was used (Schmutz et al. 
2010). 

Results 

Aphid resistance of Rasol -introduced backcross lines 

Adams and its derivatives demonstrated antibiosis to fox- 
glove. We used the term "aphid resistance" for this antibio- 
sis. The foxglove aphid resistance QTL, Rasol, localized to 
a 3.7 cM interval between Satt009 and Satt530 (Kamiya et 
al. 2008), was introduced into Toyomusume by backcross- 
ing (Fig. 1). In all TM-lines tested, the number of aphids 
per plant was significantly lower than that observed on 



SbDV tolerance of Rasol -introduced backcross lines 

Adams and its derivatives demonstrated lower incidences 
of SbDV symptoms in the field. We used the term "SbDV 
tolerance" for this partial resistance in the field. The SbDV 
tolerance of 3 independent backcross lines differing in aphid 
resistance phenotype (TM-1382, TM-1386 and TM-1394, 
highlighted in boxes in Fig. 2) was tested. Only TM-1386 
demonstrated virus tolerance equal to that of Shokukei-10, 
and TM-1382 and TM-1394 exhibited similar levels of virus 
symptom incidence as Toyomusume, the recurrent parent of 
both lines. Similar results were obtained in consecutive 
years (Fig. 3A, 3B) and when tolerance was evaluated as se- 
verity of virus symptoms (Fig. 4). 



622 



Ohnishi, Miyake, Takeuchi, Kousaka, Hiura, Kanehira, Saito, Sayama, Higashi, Ishimoto, Tanaka and Fujita 



-D 
O 



Rasol 
j (63kb) . 



B 



CO 

to 



CM 



CO 



CO 



No. of aphids 
per plant 

Exp.1 Exp.2 



Phenotype 



TM-lines 



Toyomusume [_ 
1905-179 | 



1905- 44 



1950-943 



1905-387 

1905-163 [ 

1905-602 [ 

1950-930 [ 




| 18.9** 


13.9** 


R 


| 33.4 


28.9 


s 




1 1 .9** 


R 




14.6** 


R 


I 32.4 




s 




31.6 


s 




27.3 


s 




11.8** 


R 


1 19.8** 




R 



Fig. 5. Fine mapping of Rasol to a 63-kb interval of chromosome 3. (A) Summary of the informative seven lines used to fine-map the location of 
Rasol. Line names beginning with 1905 and 1950 indicate their derivation from the mapping populations C1905 and C1950, respectively. TM- 
lines represent Toyomusume-backcrossed lines. Solid bars represent the genomic region derived from Adams, the source of aphid resistance. 
Open bars indicate the genomic region derived from Toyomusume or Toyoharuka or Shokukou-05405-1-1, the parent of the populations. Gray 
bars indicate intervals in which recombination events might occur. The scale at the top shows positions in the soybean genome sequence assembly 
(Glymal) on chromosome 3. Intervals marked with horizontal dotted lines are greatly reduced in scale. Markers used for mapping (shown at the 
bottom) are described in Table 1 . (B) Phenotype data of lines used for fine-mapping. The data represent the number of aphids per plant 14 days af- 
ter settlement. For TM-lines, the results of bulked Toyomusume-backcrossed lines are shown. Jugged phenotype by aphid numbers are drawn at 
the rightmost column. Susceptible and resistant phenotype are abbreviate to 's' and 'R' respectively. Significant differences from Toyomusume 
obtained with Dunnett's test at the 1% level are indicated by **. 



Fine mapping of Rasol 

Because Rasol solely showed the aphid resistance, fine- 
mapping of Rasol was undertaken. A total of 1245 F 2 lines 
were screened for recombination events between the mark- 
ers Satt009 and Satt530. Recombination events were con- 
finned in 93 lines, from which 42 lines were selected. The 
marker genotype for Satt641 and additional SSR markers 
located between Satt009 and Satt530 (Table 1) were deter- 
mined for the 42 lines. For genotyping, Satt641 and addi- 
tional SSR markers listed in Table 1 were used. Fifteen in- 
formative lines were selected and the degree of aphid 
antibiosis was investigated. The results for the seven most- 
informative lines is presented in Fig. 5. Rasol was located be- 
tween markers Gm03-1 1 and Gm03-12, comprising a 63 kb 
interval in the Williams 82 genome sequence. This interval 
contains three predicted genes in the Williams 82 genome 
sequence, namely Glyma03g04500 (TATA-box binding 
protein), Glyma03g04510 (calcium/calmodulin-dependent 
protein kinase-related) and Glyma03g04530 (NB-ARC do- 
main and leucine-rich repeat-containing protein, a member 
of the nucleotide-binding site-leucine-rich repeat [NBS- 
LRR] gene family) based on Phytozome v7.0. 



Discussion 

We introduced the aphid resistance QTL Rasol into a 
Toyomusume background by backcrossing. All Rasol- 
introduced lines showed the aphid resistance (Fig. 2). Hence 
we conclude Rasol represents a resistance gene(s) effective 
against foxglove aphid. However, most backcross lines 
showed a lower degree of aphid antibiosis compared with 
the donor of Rasol, Shokukei-10, or Adams. While we ob- 
served less than 7 aphids per plant 14 days after settlement in 
Shokukei-10 and Adams, we observed more than 13 aphids 
in most TM-lines (Fig. 2). These observations imply the 
presence of additional minor gene(s) that modify the degree 
of aphid antibiosis in Adams. 

In contrast to aphid resistance, introduction of Rasol by 
backcrossing was insufficient to exhibit the SbDV tolerance 
phenotype. Two backcross lines (TM-1382 and TM-1394) 
showed virus symptom incidence similar to that of the recur- 
rent parent Toyomusume, even though they harbor Rasol 
(Fig. 3). Interestingly, only one line (TM-1386) exhibited a 
SbDV-tolerant phenotype (Fig. 3). Araki etal. (1997) report- 
ed that resistance to foxglove aphids was necessary for SbDV 
tolerance in segregated progenies of Adams. Our obser- 
vations, together with those of Araki et al. (1997), suggest 
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that an additional gene(s) is required for SbDV tolerance. 
Nevertheless, because the SbDV-tolerant line TM-1386 
showed similar agricultural traits as the recurrent parent 
Toyomusume, it is useful for breeding programs in northern 
Japan. 

Rasol was fine-mapped to a 63-kb fragment containing 
three gene candidates in the Williams 82 sequence assembly 
(Glymal). One of the three gene candidates is a NB-ARC 
domain and leucine-rich repeat-containing gene (NBS- 
LRR), a typical type of so-called R-gene. In this study only 
sequence data for the soybean Williams 82 was utilized. The 
corresponding genomic sequence in Adams is unknown, 
hence we cannot discuss further. However, this NBS-LRR 
gene is a most likely candidate of Rasol because the previ- 
ously cloned aphid resistance genes Mi (Milligan et al. 1998, 
Rossi et al. 1998) and Vat (Dogimont et al. 2010), and can- 
didates of the aphid resistance genes AKR (Klingler et al. 
2005), Ragl (Kim et al. 2010a) and Rag2 (Kim et al. 
2010b), are members of the NBS-LRR gene family. 

In soybean, three resistance genes for soybean aphid 
(Aphis glycines) are reported, namely Ragl on chromosome 
7 (Kim et al. 2010a), Rag2 (Kim et al. 2010b) and Rag3 
(Zahng et al. 2010) on chromosome 13. We mapped Rasol 
on chromosome 3, therefore Rasol represents a gene inde- 
pendent from Ragl, Rag2 and Rag3. Moreover, previous 
reports demonstrate Adams does not show resistance to soy- 
bean aphid (Honda 1999), for which Ragl, Rag2 and Rag3 
show resistance. Interestingly, candidates of Ragl and Rag2 
are NBS-LRR genes, the same type of gene as one of the 
candidates of Rasol. Thus soybean might evolve novel 
NBS-LRR genes for resistance to different aphid species. 
Comparisons among Rasol, Ragl and Rag2 will help to im- 
prove our understanding of aphid resistance effected by the 
NBS-LRR gene family. 

In the case of Adams, aphid resistance and tolerance to 
the virus transmitted by the aphid are closely related. A sim- 
ilar case was reported in melon with melon-cotton aphid 
(Aphis gossypii) and Cucumber mosaic virus (CMV) 
(Boissot et al. 2010, Pitrat et al. 1980, Sarria et al. 2008). In 
this case, the causal gene Vat, which shows both aphid resis- 
tance and the resistance to aphid transmission of CMV, was 
cloned and demonstrated to be a member of the NBS-LRR 
gene family (Dogimont et al. 2010). It is interesting to deter- 
mine whether the aphid resistance manifested by NBS-LRR 
proteins will basically beside the resistance to virus trans- 
mission by the aphid or not. 
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